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Abstract 

Solving multihomogeneous systems, as a wide range of structured algebraic systems occurring fre- 
quently in practical problems, is of first importance. Experimentally, solving these systems with Grobner 
bases algorithms seems to be easier than solving homogeneous systems of the same degree. Nevertheless, 
the reasons of this behaviour are not clear. In this paper, we focus on bilinear systems (i.e. bihomoge- 
neous systems where all equations have bidegree (1,1)). Our goal is to provide a theoretical explanation 
of the aforementionned experimental behaviour and to propose new techniques to speed up the Grobner 
basis computations by using the multihomogeneous structure of those systems. The contributions are 
theoretical and practical. First, we adapt the classical F5 criterion to avoid reductions to zero which 
occur when the input is a set of bilinear polynomials. We also prove an explicit form of the Hilbert 
series of bihomogeneous ideals generated by generic bilinear polynomials and give a new upper bound 
on the degree of regularity of generic affine bilinear systems. This leads to new complexity bounds for 
solving bilinear systems. We propose also a variant of the F5 Algorithm dedicated to multihomogeneous 
systems which exploits a structural property of the Macaulay matrix which occurs on such inputs. Ex- 
perimental results show that this variant requires less time and memory than the classical homogeneous 
F5 Algorithm. 

1 Introduction 

The problem of multivariate polynomial system solving is an important topic in computer algebra since 
algebraic systems can arise from many practical applications (cryptology, robotics, real algebraic geometry, 
coding theory, signal processing, etc.). One method to solve them is based on the Grobner bases theory. 
Due to their practical importance, efficient algorithms to compute Grobner bases of algebraic systems are 
required: for instance Buchberger's Algorithm (]9], Faugere F4 IfTSl or F5 1 16|. 

In this article, we focus on the Algorithm. In particular, the F5 criterion is a tool which removes 
the so-called reductions to zero (which are useless) during the Grobner basis computation when the input 
system is a regular sequence. For instance, consider a sequence of polynomials (/i, . . . , /„). The re- 
ductions to zero come from the leading monomials in the colon ideals (/i, . . . , : /,;. Let LM(/) 
denote the ideal generated by the leading monomials of the elements of an ideal /. Then the reduc- 
tions to zero detected by the criterion are those related to LM((/i, . . . , fi-i))- For regular systems, 
LM((/i, . . . , = LM((/i, . . . , : fi). Therefore, the F5 criterion removes all useless reduc- 
tions. In practice, if a homogeneous polynomial system is chosen "at random", then it is regular 

In this paper, we consider multihomogeneous systems, which are not regular. Such systems can appear 
in cryptography JT?), in coding theory ll32l or in effective geometry (see Il35ll36l '). 
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A multihomogeneous polynomial is defined with respect to a partition of the unknowns, and is homo- 
geneous with respect to each subset of variables. The finite sequence of degrees is called the multi-degree 
of the polynomial. For instance, a bihomogeneous polynomial / of bi-degree (di, ^2) over fc[a;o, . . . , Xn^ , 
Uq, . . . , Uny] is a polynomial such that 

VA, ^, /(Axo, . . . , Xxn^ , fiyo, ^y„J = X'^^ fi'^^ f{xQ, Xn^ , 2/0, ■ • • , 2/n J- 

In general, multihomogeneous systems are not regular. Consequently, the F5 criterion does not remove all 
reductions to zero. Our goal is to understand the underlying structure of these multihomogeneous algebraic 
systems, and then use it to speed up the computation of a Grobner basis in the context of F5. In this paper, 
we focus on bihomogeneous ideals generated by polynomials of bi-degree (1,1). 



1.1 Main results 



Let fc be a field, fi, ■ ■ ■ fm G ^[a^o, ■ ■ ■ , Xn^ ,yo, ■ ■ ■ , Vny] be bilinear polynomials. We denote by Fi the 
polynomial family (/i, . . . , /,) and by li the ideal (F,). We start by describing the algorithmic results of 
the paper, obtained by exploiting the algebraic structure of bilinear systems. 

In order to understand this structure, we study properties of the jacobian matrices with respect to the 



two subsets of variables xq, . 



jacx(i^.) = 



9fi 
dxQ 



dxQ 



,Xn^ and 2/0,- 

dfi 



■ J Driy'- 



dxn 



dya 



dya 



dfi 



dy„^ 



We show that the kernels of those matrices (whose entries are linear forms) correspond to the reductions 
to zero not detected by the classical F5 criterion. In general, all elements in these kernels are vectors of 
maximal minors of the jacobian matrices (Lemma [3.1b . For instance, if = fiy = 2 and to = 4, consider 



and 



V — (minor(jac3(.(i^4), 1), — minor(jaCx(-F4), 2), minor(jaCx(-F4), 3), — minor(jac3(.(F4), 4)) 



w = (minor(jaCy(F4), 1), -minorQaCy (i^i), 2), minorQaCy (^^4), 3), -minor(jaCy (F4), 4)), 



where minor(jaCx(-F4), k) (resp. minor(jaCy (F4), fc)) denotes the determinant of the matrix obtained from 
jaCx(i^4) (resp. jaCy (F4)) by removing the fc-th column. The generic syzygies corresponding to reductions 
to zero which are not detected by the classical criterion are 

V e KerL(jaCx(^4)) and w e KeriQaCy (i^4)). 

We show (Corollarv l4.1l i that, in general, the ideal : fi is spanned by and by the maximal 
minors of jac3(.(Fi_i) (if i > Uy + 1) and jaCy (Fi_i) (if i > nx + 1). The leading monomial ideal of : 
fi describes the reductions to zero associated to fi. Thus we need results about ideals generated by maximal 
minors of matrices whose entries are linear forms in order to get a description of the syzygy module. In 
particular, we prove that, in general, grevlex Grobner bases of those ideals are linear combinations of the 
generators (Theorem 13. 2| i. Based on this result, one can compute efficiently a Grobner basis of li^i : fi 
once a Grobner basis of is known. 

This allows us to design an Algorithm (Algorithm 13. 2l l dedicated to bilinear systems, which yields an 
extension of the classical F^ criterion. This subroutine, when merged within a matricial version of the F5 
Algorithm (Algorithm l2.2l i. eliminates all reductions to zero during the computation of a Grobner basis of a 
generic bilinear system. For instance, during the computation of a grevlex Grobner basis of a system of 12 
generic biUnear equations over k[xo, . . . ,XQ,yQ, . . . ,ye], the new criterion detects 990 reductions to zero 
which are not found by the usual F5 criterion. Even if this new criterion seems to be more complicated than 
the usual F5 criterion (some precomputations have to be performed), we prove that the overcost induced 
by those precomputations is negligible compared to the cost of the whole computation. 
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Next, we introduce a notion of bi-regularity which describes the structure of generic bilinear systems. 
When the input of Algorithm l3.2l is a bi-regular system, then it returns all reductions to zero. We also give 
a complete description of the syzygy module of such systems, up to a conjecture (Conjecture 14. Il l on a 
linear algebra problem over rings. This conjecture is supported by practical experiments. We also prove 
that there are no reductions to zero with the classical F5 criterion for affine bilinear systems (Proposition 
16.11 ) which is important for practical applications. 

We describe now the main complexity results of the paper We need some results on the so-called 
Hilbert bi-series of ideals generated by bilinear systems. For bi-regular bilinear system, we give an explicit 
form of this series (Theorem l5.1l ): 

N^{tiM) = (1-^1^2)'"+ 

After this analysis, we propose a variant of the Matrix F5 Algorithm dedicated to multihomogeneous 
systems. The key idea is to decompose the Macaulay matrices into a set of smaller matrices whose row 
echelon forms can be computed independently. We provide some experimental results of an implementa- 
tion of this algorithm in Magma2 . 15. This multihomogeneous variant can be more than 20 times faster 
for bihomogeneous systems than our Magma implementation of the classical Matrix Fc, Algorithm. We 
perform a theoretical complexity analysis based on the Hilbert series in the case of bilinear systems, which 
provides an explanation of this gap. 

Finally, we establish a sharp upper bound on the degree of regularity of 0-dimensional affine bilinear 
systems (Theorem |6.1| l. Let /i, . . . , fn^+n be an affine bilinear system of k[xQ, . . . , Xn^-i, yo, ■ ■ ■ ,yn -1], 
then the maximal degree reached during the computation of a Grobner basis with respect to the grevlex or- 
dering is upper bounded by: 

dreg < min {rix + 1, Uy + 1) . 

This bound is exact in practice for generic bilinear systems and permits to derive complexity estimates for 
solving bilinear systems (Corollarv l6.1| l which can be applied to practical problems (see for instance ifTSl 
for an application to the MinRank problem). 

1.2 State of the art 

The complexity analysis that we perform by proving properties on the Hilbert bi-series of bilinear ideals 
follows a path which is similar to the one used to analyze the complexity of the F5 algorithm in the 
case of homogeneous regular sequences (see dSJ). In 1^25 J , the properties of Buchberger's Algorithm are 
investigated in the context of multi-graded rings. 

The algorithmic use of multihomogeneous structures has been investigated mostly in the framework of 
multivariate resultants (see [11. ,13 j and references therein for the most recent results) following the line 
of work initiated by ll30l . In the context of solving polynomial systems by using straight- line programs 
as data-structures, ||23l provides an alternative way to compute resultant formula for multihomogeneous 
systems. 

As we have seen in the description of the main results, the knowledge of Grobner bases of ideals 
generated by maximal minors of linear matrices play a crucial role. Theorem 13 .21 which states that such 
Grobner bases are obtained by a single row echelon form computation is a variant of the main results in 
Il38l and Q (see also the survey JS]). 

More generally, the theory of multihomogeneous elimination is investigated in (33] and f34l providing 
tools to generalize some well-known notions (e.g. Chow forms, resultant formula, heights) in the homoge- 
neous case to multihomogeneous situations. Such works are initiated in ll40l where the Hilbert bi-series of 
bihomogeneous ideals is introduced. 
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1.3 Structure of the paper 



This paper is articulated as follows. Some tools from commutative algebra are introduced. Next, we 
investigate the case of bilinear systems and propose an algorithm to remove all reductions to zero during 
the Grobner basis computation. Then we prove its correctness and explain why it is efficient for generic 
bilinear systems. To continue our study of the structure of bilinear ideals, we give the explicit form of the 
Hilbert bi- series of generic bilinear ideals. Finally, we prove a new bound on the degree of regularity of 
generic affine biUnear systems and we use it to derive new complexity bounds. Technical results and their 
proofs are postponed in Appendix. 
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2 Grobner bases: the Matrix F5 Algorithm 
2.1 Grobner bases: notations 

In this section, R denotes the ring k[xi, . . . , Xn] (where fc is a field) and for all j3 = {(3i, . . . , /3„) g N", 
denotes x^^ , ■ • • , xf^^ . Grobner bases are defined with respect to a monomial ordering (see ifTOl . page 
55, Definition 1). In this paper, we focus in particular on the so-called grevlex ordering (degree reverse 
lexicographical ordering). 

Definition 2.1. The grevlex ordering is defined by: 

{E < E A or 
E ci^i ~ E Pi '^^'^ the first coordinates 
from the right which are different satisfy ai > Pi. 

If ^ is a monomial ordering and / e i? is a polynomial, then its greatest monomial with respect to -< 
is called leading monomial and denoted by LM^(/) (or simply LM(/) when there is no ambiguity on the 
considered ordering). 

If / C i? is a polynomial ideal, its leading monomial ideal (i.e. ({LM^(/) : / G /})) is denoted by 
LM^(/) (or simply LM(/) when there is no ambiguity on the ordering) . 



5 



Definition 2.2. let I G R be an ideal, and -< be a monomial ordering. A Grobner basis of I (relatively to 
~<) is a finite subset G C I such that: (LM^(G)) = LM^(/). 

Definition 2.3. Let I d R be an ideal, -<be a monomial ordering and f R be a polynomial. Then there 
exist unique polynomials f d R and g (z I such that f ^ f + g, f is monic and none of the monomials 
appearing in f are in L M ^ (/). The polynomial f is called the normal form of f ( with respect to I and -< ), 
and is denoted NF/^^(/). 

It is well known that NF/.^(/) = if and only if / G / (see e.g. [10;]). 

Definition 2.4. Let I d R be an homogeneous ideal, ^ be a monomial ordering and D be an integer. We 
call D-Grdbner basis a finite set of polynomials G such that (G) — I and 

V/ G / with deg(/) < D, there exists g (z G such that LM-j(g) divides LM^(/). 

The following Lemma is a straightforward consequence of Dickson's Lemma |10 page 71, Theorem 

5]. 

Lemma 2.1. Let I a R be an ideal and let -<be a monomial ordering. There exists D G N such that every 
D-Grdbner basis with respect to -< is a Grobner basis of I with respect to 

2.2 The Matrix F5 Algorithm 

We use a variant of the F5 Algorithm, called Matrix F5 Algorithm, which is suitable to perform complexity 
analyses (see [4,5, 19 1). 

Given a set of generators (/i, . . . , /„) of an homogeneous polynomial ideal I C R, an integer D and 
a monomial ordering ^, the Matrix F5 Algorithm computes a D-Grobner basis of / with respect to ^. It 
performs incrementally by considering the ideals li = (/i, . . . , fi) for 1 < i < m. 

Let d G N, denote by Rii the fc- vector space of polynomials in R of degree d. As in ||T6l and ID, we use 
a definition of the row echelon form of a matrix which is slightly different from the usual definition: we call 
row echelon form the matrix obtained by applying the Gaussian elimination Algorithm without permuting 
the rows. The idea of the Matrix i^s Algorithm (see Algorithm 12. 2l below) is to calculate triangular bases 
of the vector spaces n Rd for 1 < d < D and 1 < i < m and to deduce from them a d-basis of I^+i. 
These triangular bases are obtained by computing row echelon forms of the Macaulay matrices. 

In the algorithm which follows, the columns in the matrix Aid.i correspond to the monomials of R 
of degree d and are sorted by the chosen monomial ordering -< (from the largest to the smallest). An 
homogeneous polynomial is identified with the corresponding row in the matrix. Each row has a signature 
{t, fj), where t is a monomial and 1 < j < i. The rows of the matrices are sorted as follows: a row with 
signature (^i, fj) is preceding a row with signature (^2, fk) if j < fc or (j = k and ti -< 12). 

When the row echelon form of a matrix is computed, the rows which are linear combinations of preced- 
ing rows are reduced to zero. Such computations are useless: removing these rows before computing the 
row echelon form will not modify the result but lead to significant practical improvements. The so-called 
F5 criterion (see [ 16 1) is used to detect these reductions to zero and is given below. 

Algoritlim 2.1. FsCriterion - returns a boolean 

{{t, fi) the signature of a row 
A matrix M. in row echelon form 
1: Return (t is the leading monomial of a row ofAi) 

Now, one gives a description of the Matrix Algorithm. 

Algoritlim 2.2. Matrix F5 (see P]|76l/) 

{(/i, . . . , fm) homogeneous polynomials of degree di < d2 < . . . < dm 
D an integer 
a monomial ordering -< 
Ensure: G is a D-Grdbner basis 0/ (/i , . . . , /„) for -< 
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9. 
10. 
11 
12. 
13. 
14. 
15. 
16. 
17. 
18. 
19. 
20. 
21 



for dfrom di to D do 

■Md,o "i— matrix with rows 
for ifrom \ to mdo 

Construct Add, i by adding to M.d.i~i the following rows: 
if di = d then 

add the row fi with signature {I, ft) 
end if 

if d > di then 

for all f from A4d-i,i with signature (e, fi), such that x\ is the 

greatest variable of e, add the n — A + 1 rows x\f, xx+if, . . . , x„/ with the 

signatures {xxe, fi), (xx+ie, fi), . . . , (a;„e, fi) except those which satisfy: 

Fscriterion {{xx+ke, f,),M'Z^^-l)=true 
end if 

Compute Md.i the row echelon form of Aid, i 

Add to G the polynomials corresponding to rows of A4d,i such that their 
leading monomial is different from the leading monomial of 
the row with same signature in Ai d,i 
end for 

end for 

return G 



We recall now some results mostly given by fTSl which justify the F5 criterion by relating reductions to 
zero appearing in an incremental computation of a Grobner basis of a homogeneous ideal with the syzygy 
module of the polynomial system under consideration. 

Definition 2.5. Let (/i, . . . , f„i) be polynomials of R. A syzygy is an element s ~ (si, . . . , s„j) e i?™ 
such that X] jLi fj ^3 ~ 0- ^'^^ degree of the syzygy is defined by maxj (deg(/j ) + deg(sj ) ). The set of all 
syzygies is a submodule of R"^ called the syzygy module o/(/i, . . . , /m). 

The next theorem explains how reductions to zero and syzygies are related: 

Theorem 2.1 {F5 criterion, [il6i ). 

1. Ift G LM(/i_i) then there exists a syzygy (si, . . . , Si) of (/i, . . . , fi) such that LM(si) = t. 

2. Let (t, fi) be the signature of a row of Add. m- Then the following assertions are equivalent: 

(a) the row [t, fi) is zero in the row echelon form Aid.m- 

(b) t ^ LM{Ii^i) and there exists a syzygy s = (si, . . . , s^) o/(/i, . . . , fi) such that t — LM(si). 

The rows eliminated by the criterion correspond to the trivial syzygies, i.e. the syzygies (si , . . . , Sm) 
such that VI < i < m, Si £ (/i, . . . , fi-i, fi+i, ■ ■ ■ , fm)- These particular syzygies come from the 
commutativity of R (for all I < i, j < m, fifj ~ fjfi — 0). It is well known that in the generic case, the 
syzygy module of a polynomial system is generated by the trivial syzygies. 

Definition 2.6. I\12\ page 419] Let (/i , . . . , fm) be a sequence of homogeneous polynomials and let f C R 
be the ideal (/i, . . . , fi). The following assertions are equivalent: 

L the syzygy module of (/i , . . . , fm) is generated by the trivial syzygies. 

2. for 2 < i < m, fi is not a divisor ofO in R/f^i. 

A sequence of polynomials which satisfies these conditions is called a regular sequence. 

This notion of regularity is essential since the regular sequences correspond exactly to the systems such 
that there is no reduction to zero during the computation of a Grobner basis with (see [1I6J ). Moreover, 
generic polynomial systems are regular. 
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3 Grobner bases computation for bilinear systems 



3.1 Overview 

Let F — (/i, . . . , fi) be a sequence of four bilinear polynomials in Q[xo, xi,X2, yo, 2/1,2/2], I be the ideal 
generated by F and V C he its associated algebraic variety. As above, li denotes the ideal (/i, . . . , fi), 
and we consider the grevlex ordering with xq >- ... >- Xn^ > yo >■■■>- Un ■ Since /i , . . . , /4 are 
bilinear, for all (ao, ai, 02) G and 1 < i < 4, fi{ao, ai, 02, 0, 0, 0) — 0. Hence, contains the linear 
affine subspace defined by yo = 2/1 = ?/2 = which has dimension 3. We conclude that V has dimension 
at least 3. 

Consequently, the sequence (/i, f2, fs, /a) is not regular (since the co-dimension of an ideal generated 
by a regular sequence is equal to the length of the sequence). Hence, there are reductions to zero during 
the computation of a Grobner basis with the F5 Algorithm (see 1161 ). 

When the four polynomials are chosen randomly, one remarks experimentally that these reductions cor- 
respond to the rows with signatures {xq, fi) and (?/q, fi). This experimental observation can be explained 
as follows. 

Consider the jacobian matrices 





r dh_ 

dxo 


dfi 

dxi 


dfi 1 

dx2 




dyo 


dh_ 

dyi 


dfi. 
dy2 










and jaCy(i^) = 










d.U 

- dxo 


d.U 

dxi 


dfi 

dx2 - 




. dyo 


dfi 
dyi 


dfi 
dy2 



and the vectors of variables X and Y. By Euler's formula, it is immediate that for any sequence of 
polynomials (gi , 52 , , 94), 

4 4 
(gi,---,94)-jac^(i^).X = ^gJ, and (gi, . . . , (74).jaCy (F).Y = ^ gji (1) 

i=l 1=1 

Denote by Keri(jac3j(F)) (resp. Ker^QaCy (F))) the left kernel of jac^(F) (resp. jaCy(F)). 

Therefore, if (gi, . . . , 54) belongs to Keri(jaCx(f')) (resp. Ker^QaCy (F))), then the relation ([T]l im- 
plies that ((/I, . . . , 94) belongs to the syzygy module of /. 

Given a (fc + 1, A;)-matrix M, denote by minor(M, j) the minor obtained by removing the j-th row from 
M. Consider 

V = (minor(jac^(F), 1), -minor(jaCx(F), 2), minor(jac^(F), 3), -minor(jac^(F), 4)). 

By Cramer's rule, it is straightforward to prove that v e Keri(jaCx(F)). A symmetric statement can be 
made forjaCy (F). From this observation, one deduces that minor(jaCx(-F), 4)/4 (resp. minor(jaCy (F), 4)/4) 
belongs to /s = (/l,/2,/3)■ 
We conclude that the rows with signature 

(LM(minor(jac,(F),4)),/4) and (LM(minor(jaCy(F), 4)), 74) 

are reduced to zero when performing the Matrix F5 Algorithm described in the previous section. A straight- 
forward computation shows that if F contains polynomials which are chosen randomly, then 

LM(minor(jac^(F),4)) = yjj and LM(minor(jaCy (F), 4)) = Xq. 

In this section, we generalize this approach to sequences of bilinear polynomials of arbitrary length. 
Hence, the jacobian matrices have a number of rows which is is not the number of columns incremented by 
1 . But, even in this more general setting, we exhibit a a relationship between the left kernels of the jacobian 
matrices and the syzygy module of the ideal spanned by the sequence under consideration. This allows us 
to prove a new F5 -criterion dedicated to biUnear systems. On the one hand, when plugged into the Matrix 
F5 Algorithm, this criterion detects reductions to zero which are not detected by the classical criterion. On 
the other hand, we prove that a D-Grobner basis is still computed by the Matrix F5 Algorithm when it uses 
the new criterion. 
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3.2 Jacobian matrices of bilinear systems and syzygies 

From now on, we use the following notations: 

• R = k[xo,...,Xn^,yo,...,yny]; 

• F = {fi, . . . , fm) C i?™ is a sequence of bilinear polynomials and Fi = (/i, . . . , /j) for 1 < i < 
m; 

• / is the ideal generated by F and li is the ideal generated by Ff, 

• Let M be a ^ X c matrix, with £ > c. We call maximal minors of M the determinants of the c x c 

sub-matrices of M; 

• iaCj^{Fi) andjaCy(Fi) are respectively the jacobian matrices 









r ^ .. 


dfi -1 


dxQ 


dXna, 




dyo 


dyny 






and 






dfi .. 


dfi 




dfi. . . 


dfi 


9x0 


dXnx 






dyny J 



• Given a matrix M, KerL(M) denotes the left kernel of M; 

• X is the vector of variables [xq , • . • , 

]* and Y is the vector of variables [yo, 

• (/i, . . . , fm) € k[xo, . . . , Xn^-i, Uo, • • • : is Called affine bilinear system if there exists an 
homogeneous bihnear system (/{*,..., ) e A: [xq a;„^ , 2/o 2/n„ ] such that 

fi{xo, . . .,Xnx-i,yo, ■ . .,yny-i) = fi{xo, . . . l,yo, • • • !)• 

Lemma 3.1. Let i > + 1 (resp. i > Uy + 1), and let s be a maximal minor o/jaCx(i^i-i) (resp. 
jaCy(Fj_i)J. Then there exists a vector (si, . . . , Sj_i,s) in Keri:,(jaCx(-Pi)) (resp. Keri:,(jaCy(Fi))J. 

Proof. The proof is done when considering s as a maximal minor of jaCx(-Fi_i) with i > + 1. The case 
where s is a maximal minor of jaCy with i > rij^ + 1 is proved similarly. 

Note that jaCx(-Fi-i) is a matrix with i — 1 rows and + 1 columns and i — 1 > + 1. Denote by 
(ii) ■ • ■ 7 the rows deleted from jac3(.(Fj_i) to construct its submatrix J whose determinant is s. 

Consider now the i x [i — Ux — 2)-matrix T such that its {£, k) entry is 1 if and only if ^ = jk else it is 
0. N denotes the following i x (i — 1) matrix: 

N = [ jacjF,) I T ] . 

A straightforward use of Cramer's rule shows that 

(minor(N,l),-minor(N,2),...,(-l)'+^minor(N,i)) e KerL(N). 

Remark that this imphes 

(minor(N, 1), -minor(N, 2), . . . , (-l)*+iminor(N, i)) G KerL(jac^(Fi)). 

A routine computation of minor(N, z) by going across the last columns of N shows that minor(N, i) = ±s 

□ 

Theorem 3.1. Let i > + 1 (resp. i > Uy + 1) and let she a linear combination of maximal minors of 
jaCx(-Fj-i) (resp. jaCy(i=i_i);. Then s e li-i : fi. 
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Proof. By assumption, s = J^i'^e^e where each Si is a maximal minor of iac^{Fi^i). According to 
Lemma im for each minor Se there exists {sf^ , • ■ • , such that 

Thus, by summation over £, one obtains 

ais['\ . . . , ^ a41,,s) e KerL(jac,(F,)). (2) 

I e 

Moreover, by Euler's formula 

By the relation ©, s /, + Y.'i2x (E^ 

OfS,- ) /j = 0, which implies that s e : /i. □ 

Corollary 3.1. Let i > n^; + 1 f resp. i > rtj, + 1 j, A/x*' f reip. A/y*' j he the ideal generated by the maximal 
minors of ]ac^{Fi) (resp. jaCy(F^)j. Then M^'^^^ C li-i : ft (resp. M^^^^^ C h-i : ft). 

Proof. By Theorem 13. II all minors of jaCx(^i-i) (resp. jaCy(Fi_i)) are elements of : fi. Thus, 

li-i : fi contains a set of generators of ^"^ (resp. My ^^). Since : fi is an ideal, our assertion 
follows. □ 
The above result implies that for all g e (resp. g e Aly ^'), the rows of signature (LM(g), fi) 

are reduced to zero during the Matrix F5 Algorithm. In order to remove these rows, it is crucial to compute 
a Grobner basis of the ideals A/x and My . These ideals are generated by the maximal minors of 
matrices whose entries are linear forms. The goal of the following section is to understand the structure of 
such ideals and how Grobner bases can be efficiently computed in that case. 



3.3 Grobner bases and maximal minors of matrices with linear entries 

Let ^ be the set of homogeneous linear forms in the ring i?x — k[xo, . . . , Xn^ ] , ^ be the grevlex ordering 
on i?x (with xq y ■ ■ ■ >■ Xn^) and Mat_5f (p, q) be the set of p x g matrices with entries in ^ with p > q 
and Tlx > p — q. Note that Mat^(p, q) is a fc-vector space of finite dimension. 

Given M e Mat^(p. q), we denote by MaxMinors(M) the set of maximal minors of M. We denote by 
Macaulay^ (MaxMinors(M), q) the Macaulay matrix in degree q associated to MaxMinors(M) and to the 
ordering -< (each row represents a polynomial of MaxMinors(9n) and the columns represent the monomials 
of degree q of k[xo, . . . , XnJ sorted by -< from the largest to the smallest). 

The main result of this paragraph lies in the following theorem: it states that, in general, a Grobner 
basis of (MaxMinors(M)) is a linear combination of the generators. 

Theorem 3.2. There exists a nonempty Zariski-open set O in Mat^{p, q) such that for all W\ E O, a 
grevlex Grobner basis of (MaxMinors(M)) with respect to -< is obtained by computing the row echelon 
form o/Macaulay^(MaxMinors(M), q). 

This theorem is related with a result from Sturmfels, Bernstein and Zelevinsky (1993), which states that 
the ideal generated by the maximal minors of a matrix whose entries are variables is a universal Grobner 
Basis. We tried without success to use this result in order to prove Theorem 13. 21 Therefore, we propose 
an ad-hoc proof, which is based on the following Lemmas whom proofs are postponed at the end of the 
paragraph. 

Lemma 3.2. Let MonomialSp_g((7) be the set of monomials of degree q in k[xQ^ . . . , Xp_q]. There exists a 
Zariski-open subset O' o/Mat_5f (p, q) such that for all M G O' 

(MonomialSp_q(g)) C LM((MaxMinors(M))) 
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Lemma 3.3. Let MonomialSp_q(g) be the set of monomials of degree q in k[xo, ■ ■ ■ , Xp^q]. There exists a 
Zariski-open subset O" q/Mat^(p, q) such that for all M G O" 

LM((MaxMi nors(M))) C (MonomialSp_q((7)) 

Lemma 3.4. The Zariski-open set O' D O" C Mat^(p, q) is nonempty. 

Proof of Theorem \3.2\ From Lemmas |3 . 21 13 . 3 1 and 13 .41 O = O' n O" is a nonempty Zariski open set. Now 
let M be a matrix in O C Mat_^(p. q). 

(MonomialSp_q(g)) = LM((MaxMinors(M))). 

Thus all polynomials in a minimal Grobner basis of (MaxMinors(M)) have degree q and then can be 
obtained by computing the row echelon form of Macaulay,^ (MaxMinors(M), q). □ 

We prove now Lemmas 13.21 13. 3l and l3. 41 
Proof of Lemma \T2\ Let 9Jl be the {p, q)-mat[ix whose (i, j)-entry is a generic homogeneous linear form 



Denote by 



a ^ {a^k'^\0 <k<n^, I < i < p, 1 < j < q} 



and given a set 

a = {a[^''^' e k,0 < k <n^, 1 < i < p, 1 < j < q} 
consider the specialization map (p^ : 'OJl i-^ 3Jla G Mat^{p, q) such that the (i,j)-entry of 97ta is 
J2k=o ^fc '"'^^fc ^ fc[a;o: • • ■ : XnJ- Wc prove below that there exists a polynomial g G k[a] such that, if 
5(a) 7^ then 

(MonomialSp_,(g)) C LM((MaxMinors((/?a(9^))))- 
Consider the Macaulay matrix Macaulay^(MaxMinors(EOt), q). 

Remark that the number of monomials in MonomialSp_g((7) equals the number of maximal minors of 
9JT. Moreover, by construction of Macaulay^(MaxMinors(9Jt), q) and by definition of -< (see Definition 
12. It , the first columns of Macaulay^ (MaxMinors(S[Jl), q) contain the coefficients of the monomials in 
MonomialSp_g(g) of the polynomials in MaxMinors(9Jt). 

Saying that (MonomialSp_g(g)) C LM((MaxMinors(9Jt))) is equivalent to saying that the determinant 
of the square submatrix of Macaulay^ (MaxMinors(2Jl), q) containing its first (p columns is non-zero. Let 
g G k[a] be this determinant. 

The inequation 5 7^ defines a Zariski-open set O' such that for all a G O' 



(MonomialSp_,((7)) C LM((MaxMinors((y5a(9K)))). 

In the following denotes the canonical inclusion morphism from k[xo, . . . , Xn^ 
where k' is the field of fractions fc(a;p_g+i , . . . , ). 



□ 



to k' [xq , 



For {vi , 



-0V : 



), denotes the specialization morphism: 



fc[a;o, ■ • ■ , Xn^\ ^ fc[iEo, • • • , Xp—q\ 

f {XO: • • • : Xnj.^ I ^ /(xq , • . • , Xp_g , 1^1 , . . . , Z^^^ _p-f g) 



Lemma 3.5. There exists a Zariski open set O'", such that i/a G O'", then the ideal (MaxMinors(?/; • 
(y5a(S0^))) w radical and its degree is {^^1)- 

Proof. There exists an affine bilinear system /i, . . . , /p G fc'(a)[a;o, . . . , Xp-q, yo, . . . , yq-2], such that: 



( \ 



yq-2 

V 1 / 
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Let / denote the ideal (/i, . . . , fp). According to Lemma IB. 31 (in Appendix), there exists a polyno- 
mial hi G k[a], such that if hi{a) ^ 0, then ^(MaxMinors(7/; o (p^{m))) (<(5a(/i), ■ ■ ■ , fMp)) 

h [xq , . . . , Xp—q] . 

One remarks that there also exists a polynomial h2 G k[a] such that if /12(a) 0, then (fiail) is 
0-dimensional (since /i , . . . , /p is a generic affine bilinear system with p equations and p variables, see 
Proposition IA.3b . From Lemma IB721 (in Appendix), there exists a polynomial such that if /13(a) 7^ 
0, then (pa.{I) is radical. From now on, we suppose that /ii(a)/i2(a)/i3(a) 7^ 0. If {wq, . . . ,Wp-q) € 
Far((MaxMinors('i/' o (pa(2t)))) (where Var denotes the variety), then the set of points in Var{(pg^{I)) 
whose projection is {wq, . . . , Wp^q) can be obtained by solving an affine linear system. The set of solutions 
of this system is nonempty and finite (since (pa{I) is 0-dimensional), thus it contains a unique element. So 
there is a bijection between Var{(pa{I)) and l/ar((MaxMinors(^/; o ipi^(9Jl)))). Since fs,{I) is radical, 

deg(v3a(/)) = deg(\/(MaxMinors(V' o Lp^i^)))). 
From Corollarv lB.il this degree is ( According to Lemma [J!2l 



deg(^(MaxMinors('(/' o ^a(9K)))) < dcg((MaxMinors(V' o <^a(9K)))) 

< dcg((MonomialSp_g((z))) = (/i). 

Therefore, 

deg(v/(MaxMinors(V' o <^a(9K)))) = deg((MaxMinors('0 o Lpa{m)))) 

and thus 

V(MaxMinors('0 o <Pa(SW))) = (MaxMinors(?A o ip^{m))). 

Furthermore, the inequation /ii (a) /12(a) /13(a) 7^ defines the wanted Zariski open set. □ 

Proof of Lemma U3\ Consider the Zariski open set O" = O' n O'" (where O' is defined in Lemma [l!2l and 
O'" is defined in Lemma [33] l and let a be taken in O" . According to Lemma [J!2l 

MonomialSp_,(gf) C LM((MaxMinors(-0 o (/Pa(9Jl)))). 

A basis of fc'[xo, . . . , a;p_^]/(MonomialSp_^(g)) is given by the set of all monomials of degree less than q. 
Therefore, the dimension of fc'[xo, . . . , Xp^q]/ (MonomialSp_g(g)) (as a /c'-vector space) is {^^-^- Thus, 
from Lemma [331 



deg((MaxMinors(V'o^a(a'T)))) = 




= deg((MonomialSp_g((7))). 



Therefore, all polynomials in (MaxMinors('(/; o (pa(2t))) have degree at least q. 

Now let g 7^ be a polynomial in (MaxMinors((pa(2t)))- Then there exists v = (wi, . . . , Vn^~p+q) 
such that the specialized polynomial verifies ipv{g) 7^ and such that deg((MaxMinors(?/'v ° ¥'a(9^)))) = 
(^^j). Thus ipv{g) is a polynomial of degree at least q in k[xo, . . . ,Xp-q]. Now suppose by con- 
tradiction that LM(g) ^ (MonomialSp_g(g)). Since deg{ipv{g)) > q, there exists a monomial m in 
g such that m G (MonomialSp_q(q)). Thus consider gi — g — Am + ANF(m). One remarks that 
LM{g) = LM(,gi) ^ (MonomialSp_g((7)). Since gi G (MaxMinors((^a(9K))), by a similar argument 
there also exists a monomial mi G (MonomialSp_g(g)) in gi. By induction construct the sequence 
gi = gi-i — Ai_imi_i + Ai_iNF(mi_i). This sequence is infinite and strictly decreasing (for the in- 
duced partial ordering on polynomials: hi -< /i2 if LM(/ii) -< LM(/i2) or if LM(/ii) — LM(/i2) and 
hi — LM(/ii) -< /i2 — LM(/i2)). But, when -< is the grevlex ordering, there does not exist such an infinite 
and strictly decreasing sequence. 

Therefore LM (5) G (MonomialSp_g((7)), which concludes the proof. □ 
Proof of Lemma \3~4\ In order to prove that the Zariski open set O' D O" is nonempty, we exhibit an 
explicit element. Consider the matrix M of Mat_5f(p,(7) whose (i, j)-entry is a;i+j_2ifO < i+j — 2 <p—q 
and i > j, else it is 0. 
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/ xo . . . \ 

Xi Xq ' ■ 

Xp—q—i 

V ... Xp-g J 

Remark that MaxMinors(M)) C k[xQ, . . . ,Xp^q]. Since (MonomialSp_q(g)) is a zero-dimensional 
ideal in k[xo, . . . ,Xp^q], the fact that LM(MaxMinors(M)) = MonomialSp_q((7) implies the equality of 
the monomial ideals LM((MaxMinors(M))) = (MonomialSp_q(q)). Thus, we prove in the sequel that 
LM(MaxMinors(M)) = MonomialSp_g((7). 

A first observation is that the cardinality of MaxMinors(M) equals the cardinality of MonomialSp_q((7). 
Let m be a maximal minor of M. Thus m is the determinant ofaqxq submatrix M' obtained by removing 
p — q rows from M. Let ii, . . . , ip-q be the indices of these rows (with ii < . . . < ip-^q). Denote by * 
the product coefficient by coefficient of two matrices (i.e. the Hadamard product) and let &q be the set of 
qy. q permutation matrices. Thus m — X^cree (— l)^^"''^-' det(cr * M'). 

Since for all a G &q, det(cr * M') is a monomial, there exists cr" £ &q such that LM(to) = ± det((T° * 
M'). 

We prove now that (t° = id. Suppose by contradiction that cj" 7^ id. In the sequel, we denote by 

• M'[i,j] the (i, j)-enti-y of M'. 

• Gi the q X 1 unit vector whose i-th coordinate is 1 and all its other coordinates are 0; 

• aj is the integer i such that cr"ej = e,;. 

Since, by assumption, cr*' 7^ id, there exists 1 < i < j < q such that ct" > af. Because of the structure of 
M, we know that for the grevlex ordering xq )~ ■ ■ ■ >- x„^, 

M'[i,aO]M'[j,a,"]^M'[z,a°]M'b-,aO]. 

Let a' be defined by 

{a^if k ^ i and k ^ j 
a'^ifk^t 
alxfk^j 

Thendet(cr'*M') >■ det(cr° * M') and by induction det(id * M') >- det(cr° ★ M'). This also proves that the 
coefficient of det (id *M') in MaxMinors(M) is 1 and contradicts the fact that LM(m) = ± det((T° ★ M'). 
This proved that LM(m) — \ det(id * M')|. Now one can remark that 

det(id * M') = xl'-^x^^-''~\'i~''~^ . . . xlZi[-''\ 

If mi, m2 are distinct elements in MaxMinors(M), then LM(mi) 7^ LM(to2)- For all to in MaxMinors(M), 
LM(to) e MonomialSp_g((7), and MaxMinors(M) has the same cardinality as MonomialSp_q((7). There- 
fore, one can deduce that LM(MaxMinors(M)) = MonomialSp_g((7). □ 

3.4 An extension of the F5 criterion for bilinear systems 

We can now present the main algorithm of this section. Given a sequence of homogeneous bilinear forms 
F — (/i, . . . , /,„) C R generating an ideal I C R, ^ the grevlex monomial ordering on R with xo >- 
■ ■ ■ Xn^ >- 2/0 >- • • • 2/ri ^ it returns a set of pairs (g, fi) such that g G li-i : fi and g ^ (for 
i > min(na; + l,ny + 1)). Following Theorem 13 . 1 1 and |3 .21 this is done by considering the matrices 
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Require: 



}ac^{Fi) (resp. jaCy(Fi)) for i > n.j; + 1 (resp. i > Uy + 1) and performing a row echelon form on 
Macaulay^(MaxMinors(jaCx(Fi)), ria; + 1) (resp. Macaulay^(MaxMinors(jaCy(Fi)), rij, + 1)). 

First we describe the subroutine Reduce (Algorithm l3.1b which reduces a set of homogeneous polyno- 
mials of the same degree: 

Algorithm 3.1. Reduce 

Require: {S, q) where S is a set of homogeneous polynomials of degree q. 
Ensure: T is a reduced set of homogeneous polynomials of degree q. 

1: M Macaulay(5, q). 

2: M ^ RowEchelonForm(M). 

3: Return T the set of polynomials corresponding to the rows o/M. 

The main algorithm uses this subroutine in order to compute a row echelon form of the matrix 
Macaulay^(MaxMinors(jaCx(-Pi))i "a: + 1) (resp. Macaulay^(MaxMinors(jaCy(Fi)), rij, + 1)): 

Algorithm 3.2. BLcriterion 

m bilinear polynomials /i, . . . , such that m < + Hy. 
< a monomial ordering over k[xQ, . . . , a;„^ ,yo, ■ ■ ■ , yuy] 
Ensure: V a set of pairs (h, fi) such that h G li-i : fi- 
V ^% 

for ifrom 2 to m do 
if i > Uy + 1 then 

T ^ Reduce(MaxMinors(jaCy(i^i_i)),r7,y + 1). 
for h in T do 

V^V\J{{hJ,)} 
end for 
end if 

9: if i > Tlx + 1 then 

10: T' Reduce(MaxMinors(jaCx(-Fi-i)), n,j: + 1). 

ll: for h in T' do 

12: V^V\j{{hJ,)} 

13: end for 

14: end if 

15: end for 

16: Return V 



The following Proposition explains how the output of Algorithm 13.21 is related to reductions to zero 
occurring during the Matrix Fc, Algorithm. 

Proposition 3.1 (Extended Fc, criterion for bilinear systems). Let fi, . . . , fm be bilinear polynomials and 
-< be a monomial ordering. Let {t, fi) be the signature of a row during the Matrix F^ Algorithm and let V 
be the output of Algorithm BLCRITERION. Then if there exists {h, fi) in V such that LM{h) = t, then the 
row with signature (t, fi) will be reduced to zero. 



Proof. According to Theorem l3.ll hfi e h-i. Therefore 

i-l 

tfr = {h~t)f,+Y,gjfj. 

i=i 

This implies that the row with signature (i, fi) is a linear combination of preceding rows in the matrix 
Macaulay(i^i, deg{tfi)). Hence this row will be reduced to zero. □ 

Now we can merge this extended criterion with the Matrix F5 Algorithm. To do so, we denote by V the 
output of BLcriterion {V has to be computed at the beginning of Matrix F^ Algorithm), and we replace 
in Algorithmic] the i^sCRlTERlON by the following BlLlNi^s CRITERION: 
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4 Fs without reduction to zero for generic bilinear systems 

4.1 Main results 

The goal of this part of the paper is to show that Algorithm 13.21 finds all reductions to zero for generic 
bilinear systems. In order to describe the structure of ideals generated by generic bilinear systems, we 
define a notion of bi-regularity (Definition l4.1l ). For bi-regular systems, we give a complete description of 
the syzygy module (Proposition |4.2| and Corollary 14. It . Finally, we show that, for such systems. Algorithm 
|3.2| finds all reductions to zero and that generic bilinear systems are bi-regular (Theorem |4.1| l, assuming a 
conjecture about the kernel of generic matrices whose entries are linear forms (Coniecture |4.1| ). 

4.2 Kernel of matrices whose entries are linear forms 

Consider an monomial ordering -< such that its restriction to k[xo, . . . , a;„^] (resp. k[yo, . . . , yuy]) is the 
grevlex ordering (for instance the usual grevlex ordering with xq xi . . . ^ . . . ^ Vny)- 

Let £,c,nx be integers such that c < £ < + c — \. Let M be the set of matrices £ x c where 
coefficients are linear forms of fc[a;o, . . . , XnJ\- Let T be the set of £ x — c — 1) matrices T such that: 

• each column of T has exactly one 1 and the rest of the coefficients are 0. 

• each row of T has at most one 1 and all the other coefficients are 0. 

• (T[ii, ji] = T[i2, j2] = 1 and ii < 12) ji < j2 

If T e T and M G 7M, we denote by Mj the £ x — 1) matrix obtained by adding to M the columns of 
T. According to the proof of Lemma [3T1 some elements of the left kernel of a matrix M can be expressed 
as vectors of maximal minors: 



Actually, we observed experimentally that kernels of random matrices M G are generated by those 
vectors of minors. This leads to the formulation of the following conjecture: 

Conjecture 4.1. The set of matrices M G such that 



VTgT, 



minor(MT,l) \ 
-minor(MT,2) 



G KerL(M). 



V(-i)' 



minor(MT, m) J 




I 



minor(MT, 1) 
-minor(MT, 2) 



> 




minor(MT, rn) ) 



contains a nonempty Zariski open subset ofM. 
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4.3 Structure of generic bilinear systems 

With the following definition, we try to give an analog of regular sequences for bilinear systems. This 
definition is closely related to the generic behaviour of Algorithm 1 3. 2 1 

Remark 4.1. In the following, Monom ia Is^ (d) fresp. Mor\om\a\s^{d)) denotes the set of monomials of 
degree d in k\xo^ . . . , Xn] (resp. fc[?/o, ■ • ■ , 2/n]j- If n < 0, we use the convention Monomials^((i) = 
Monomials^(d) = 0. 

Definition 4.1. Let m < nx + riy and /i , . . . , be bilinear polynomials of R. We say that the polynomial 
sequence (/i, . . . , /,„) a bi-regular sequence //rn — 1 orif{fi, . . . , fm-i) is a bi- regular sequence and 

LM(/„_i : /„) = (Monomials^_„^_2(% + 1)> 
+ (Monomialsf„_„^_2("x + 1)> 
+ LM(/™_i) 

In the following, we use the notations: 

• B£{nx, riy) the fc-vector space of bilinear polynomials \n K[xa, . . . , Xn^ ,yo, ■ ■ ■ , Uuy]', 

• X C k[xo, . . . , Xn^ , 2/0, • • ■ , J/n^] (rcsp. Y) is the ideal ( ) (resp. {yo,...,y„jy, 

• An ideal is called bihomogeneous if there exists a set of bihomogeneous generators. In particular, 
ideals spanned by bilinear polynomials are bihomogeneous. 

• Ji denotes the saturated ideal : {X nY)°°; 

• Given a polynomial sequence (/i, . . . , /„,,), we denote by Syztriv the module of trivial syzygies, i.e. 
the set of all syzygies (si, . . . , s„i) such that 

• A primary ideal P <Z Ris called admissible if {xq, ■ ■ ■ , XnJ) 't- ^/P ™d (yo, ■ • ■ , yny) ^ ^/P\ 

• Let E' be a fc-vector space such that dim(i?) < oo. We say that a property V is generic if it is satisfied 
on a nonempty open subset of E (for the Zariski topology), i.e. 3/i e fc[oi, . . . , adiin(£;)], h ^ 0, 
such that 

V is does not hold on (oi, . . . ,adim(£;)) =^ H^^i, ■ ■ • ,adiin(B)) = 0. 
Without loss of generality, we suppose in the sequel that < Uy. 

Lemma 4.1. Let be an ideal spanned by m generic bilinear equations /i , . . . , /„ and Im — PevP 
be a minimal primary decomposition. Let Po ^ V be one of its primary non-admissible components. If 
m < Tlx + 1 (resp. m < Uy -\- 1), then X <f_ VPq (resp. Y <f_ VPq). 

Proof. Suppose that m < Ux + 1- Consider the field k' — k{yo, . . . , ?/„^) and the canonical inclusion 

ijj : R^ k'[xQ, . . . , x„J. 

ip{Im) is an ideal of fc' [a;o , . . . , ] spanned by m polynomials of k'[xQ, . . . , Xn^ ] ■ Generically, the system 
(■0(/i), . . . , ipifm)) is a regular sequence of k'[xQ, . . . , Xn^]- Thus there exists an polynomial f E X 
(homogeneous in the XiS) such that ■ip{f) is not a divisor of in k'[xo, . . . , Xn^]/4'{Irn)- This means that 
4>{Im) '■ ^(/) = 4'{Im)- Suppose the assertion of Lemma 14. II is false. Then X C \/Pq and hence, / € 
\/Pq. Therefore there exists g e fc[yo, • ■ • , ?/n„] such that, in R, gf e V7m (take g in (npg-py^p^j-v/P) \ 
{^/Po} which is nonempty). Thus ip{f) £ \/ V'(-^m) (since '(/'(ff) is invertible in fc'), which is impossible 
since : = ^ijm)- □ 

Lemma 4.2. * If "m < Ux there exists a nonempty Zariski-open set O C BCK{nx,ny)™ such that 
(/i , . . . , fm) C O implies that has co-dimension m and all the components of a minimal primary 
decomposition of Im are admissible; 
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• if Hx + i < rn, then there exists a nonempty Zariski-open set O C BCKinx,ny)"^ such that 
(/i, . . . , fm) C O implies that X is a prime associated to \J Im', 

• ;/ riy + 1 < m, then there exists a nonempty Zariski-open set O C BCK{nx,ny)™ such that 
(/i, ■ • ■ , /m) C O implies that Y is a prime associated to vXn- 



Proof. 

• If TO < 71a;, then by Lemma l4n J,n = Im- Then according to Theorem lA.il there exists a nonempty 
Zariski-open set O C BCxinx, ny)"^ such that (/i, . . . , /,„) C O impHes that (/i, . . . , /,„) is a 
regular sequence. Therefore, /,„ has co-dimension m and all the components of a minimal primary 
decomposition of /,„ are admissible. 

• If Tlx + I < m, then according to Proposition I A. 31 J„i = {I„i : Y°°) : X°° is equidimensional of 
co-dimension to. Let Vx be the set {(0, . . . , 0, ao, . . . , a„ )\ai £ k}. Since Vx C Var{Im : Y°°) 
and codim(T4) = Ux + l, it can be deduced that Vx 'f. Var{J,n) and Var{Im '■ Y°°) = Var(Jm) U 
Vx. This means that ^/Im ■ Y°° = \/Jm H X and yOm 't- Thus X is a prime associated to 
\/Ira ■ Y°°. Since Y is not a subset of X, X is also a prime ideal associated to ^/T^. 

• Similar proof in the case n^, + 1 < to. 

□ 

Lemma 4.3. Suppose that the local ring Rx /Ix ( resp. Ry /ly ) is regular and that X ( resp. Y) is a prime 
ideal associated to \fl and let Q be an isolated primary component of a minimal primary decomposition 
of I containing X (resp. Y). Then Q = X (resp. Q = Y). 

Proof. By assumption, X is a prime ideal associated to \/7- Then, there exists an isolated primary com- 
ponent of a minimal primary decomposition of / which contains a power of X and does not meet R \ X. 
This proves that Ix does not contain a unit in Rx ■ 

By assumption Rx/Ix is regular and local, then Rx/Ix is an integral ring (see e.g. lfT2l Corollary 
10.14]) which implies that Ix is prime and does not contain a unit in Rx- 

Let / — Qin-'-nQsbea minimal primary decomposition of /. In the sequel, Qi^ denotes the 
localization of Qi by X. Suppose first that there exists 1 < i < s such that Ix — Qix with Qi non- 
admissible which does not meet the multiplicatively closed part R \ X . Then Qi^ is obviously prime 
which implies that Qi itself is prime [3, Proposition 3.11 (iv)]. Our claim follows. 

It remains to prove that Ix = Qix for some I < i < s. Suppose that the Qj's are numbered such that 
meets the multiplicatively closed set for r+l <j< s but not Qi, . . . , Q^- Ix ~ Qix^' ' '^Qrx 

and it is a minimal primary decomposition |3 , Proposition 4.9]. Hence, since Ix is prime, r — 1 and Qi is 

the isolated minimal primary component containing X. 

Proving that Q = Y in the case where Ry / ly is regular and that F is a prime associated to %// is done 

in the same way. □ 

Proposition 4.1. Let k be a field of characteristic 0. There exists a nonempty Zariski-open set O C 
BC{nx, ny)™ such that for all (/i, . . . , /,„) C O the non-admissible components of a minimal primary 
decomposition of {fi, . . . , fm) are either X or Y. 



Proof. Suppose that Ux + 1 < m. Then, from Lemma 1431 there exists a nonempty Zariski-open set Oi 
such that X is an associated prime to \/T- Note also that this implies that Ix has co-dimension Ux + 1. 
Thus, from Lemma |431 it is sufficient to prove that there exists a nonempty Zariski-open set O2 such that 
for all (/i, . . . , fm) & Old O2, Rx/Ix is a regular local ring. 

From the Jacobian Criterion (see e.g. |12], Theorem 16.19), the local ring Rx/Ix is regular if and 
only if jac(/i, . . . , fm) taken modulo X has co-dimension Hx + 1. Since the generators of / are bilinear, 
the latter condition is equivalent to saying that the matrix 



Jx 



dfi ... _dfi_ 

dxQ dxn^ 



dfm. . . . d.U 
dxQ dxn^ 
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has rank Ux + I- We prove below that there exists a nonempty Zariski-open set O3 such that for all 
(/i, . . . , fm) e O3, Jx has rank + 1. 

Let Ci, . . . , c„i be vectors of coordinates of B£{nx, ny)"^, DJl be the vector of all bilinear monomi- 
als in R with respect to the partition [a;o, • ■ • , Xn^], [uo, ■ ■ ■ , Uuy] and ^ be the field of rational fractions 
A:(ci, . . . , Cm)- Consider the polynomials fi = DJl.cJ for 1 < i < m and the Zariski-open set O3 in 
B£{nx, ny)™ defined by the non-vanishing of all the coefficients of the maximal minors of the matrix 



dxo 



dxo 



J3fi 



dXn 



It is obvious that (/i, . . . , /„) e O3 implies that Jx has rank rix + l; our claim follows. 

In the case where Uy < m. The proof follows the same pattern using Lemmas 14.21 and 14.31 and the 
Jacobian criterion. The only difference is that one has to prove that there exists a nonempty Zariski-open 
set O4 such that for all (/i ,...,/„) e O4 the matrix 



Jy 





dfi 1 


dya 








dyo 





has rank Uy + 1, which is done as above. □ 

Remark 4.2. The proof of Proposition \4.1\ relies on the use of the Jacobian Criterion. From fJ2, Theorem 
16.19], it remains valid if the characteristic ofk is large enough so that the residue class field of X (resp. 
Y) is separable. 

The two following propositions explain why the rows reduced to zero in the generic case during the 
Algorithm have a signature (t, fi) such that t e k[xo, . . . , XnJ\ or t S fc[?/o, ■ • • , Tjuy]- 

Proposition 4.2. Let m be an integer such that m < Ux + Uy. Let L be the set of bilinear systems 
with m polynomials (L (Z R"^). Then the set of bilinear systems /i, . . . , such that Syz = {{Syz fl 
k[xo, ■ . ■ , a^n^]™) U {Syz D k[yQ, . . . , 2/nj,]™) U Syztriv) is a nonempty Zariski-open subset of L. 

Proof. Let s = (si, . . . , Sm) be a syzygy. Thus, Sm is in /m-i : fm- We can suppose without loss of 
generality that the Si are bihomogeneous of same bi-degree (Proposition lA.ll ). According to Theorem 
lA.ll there exists a nonempty Zariski open set Oi C BC{nx,ny)"^, such that if (/i, . . . , fm) S Oi, 
then /,„ is not a divisor of in R/Jm~i- We can deduce from this observation that Sm € Jm-i- So 
Sm € Im-1 or there exists P a non-admissible primary component of Im-i such that s,„ ^ P. Assume 
that Sm ^ Im-i- From Proposition l4.1l there exists a nonempty Zariski open set O2 C B£{nx,ny)"\ such 
that if e O2, then ( ) = P (or (j/o, • • • , Vjiy = P)- This means that, generically, 

Sm e k[yo, y„J (or s,„ e fc[xo, . . . , XnJ). 

Finally, we see that, if (/i, . . . , /,„) e Oi n O2, then s,„ e U k[yo, ?/„J U k[xo, . . . , x„J. 

Since the syzygy module of a bihomogeneous system is generated by bihomogeneous syzygies, it can be 
deduced that Syz = {{Syz n k[xo, x„J™) U {Syz n k[yo, . . . , y„J™) U Syztriv)- □ 

Proposition 4.3. Let V be the output of Algorithm BLCRITERION and let {h, fi) be an element ofV. Then 

• ifh e k[xo, . . . , XnJ, then Vj, yjh € 

• ifh e /c[yo, . . . , y„J, then Vj, Xjh e Ii-i. 

Proof. Suppose that h E k[xo, . . . , XnJ is a maximal minor of jaCy (the proof is similar if h E 

k[yo, ■ ■ ■ ,yn ])■ Consider the matrix ]ac^{Fi^i) as defined in Algorithm 13.21 Then there exists an {i — 
1) X {i — 1) extension Mt of jaCx(-Fi_i) such that det(AfT) — h (similarly to the proof of Lemma [TTT l. 
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Let < j < rij, be an integer. Consider the polynomials hi, ... , /i^-i, where hk is the determinant of the 
{i — 2) X {i — 2) matrix obtained by removing the (j + l)th column and the fcth row from Mt- 
Then we can remark that 

{hi ~h2 ... (-l)'-i/i,_2 • Mt = (0 ... (-IpdetCMr) ... O) 

where the only non-zero component is in the [j + l)th column. Keeping only the Uy + 1 first columns of 
Mt, we obtain 

{hi -h2 ... (-l)"«/i„„+i) •jac,(F,_i) = (0 ... {-lydetiAr) ... O) 

/i 

Since jac^(Fi_i) • ■ I = I I , the following equality holds 





{hi ~h2 ... {-ly^y-^hn^ (~l)"«;i„^+i) . : \ ^y^det{MT)^yjh. 



This implies that yjh £ li^i. □ 

Corollary 4.1. Let m be an integer such that m < tlx + Uy and let fi, . . . , Jm be bilinear polynomials. 
Let V be the output of Algorithm BLCRITERION. Assume that 

(/.,„_! : /™) n fc[xo, ... ,a;„^] = ({/i e fc[xo, ... ,a;„^] : {hJyn)^V}). 

(^i-i : /m) n fc[?/o, • ■ ■ ,ynj = ({/i e fc[2/o, ■ • • ,ynj : {h,!„i)^V}). 

LetGx (respGy)beaGrdbnerbasisof{Im-i ■ fm)C^k[xo, . . . ,XnJ\(resp. (/,„_i : /m)nfc[2/o, • • ■ , J/n^lJ 
and let Gm-i bea Grdbner basis of Im-i- IfSyz — {{Syzr]k[xo, . . . , a;„^]™)U(S'yznfc[?/o, ■ • • i2/r!„]'")U 
Syztriv), then Gx U Gy U G„i-i is a Grdbner basis of Im-i ■ fm- 

Proof Let / e Im-i ■ fm be a polynomial. Thus there exist si, . . . , s„j_i such that (si, . . . , Sm-i, /) G 
Syz. Since Im-i and fm are bihomogeneous, we can suppose without loss of generality that / is bihomo- 
geneous (Proposition lA. 11 ). Let (di, ^2) denote its bi-degree. 

• If da = (resp. di = 0), then / e {G^) (resp. / e (Gy)). 



Let Gx = {5i'^''}i<t<card(G,) and Gy = {gi^'}i<i<card(G„)- If rfi 7^ and ^2 ^ then, since 
Syz = {{Syz n k[xo, x„J") U {Syz n k[yo, . . . , 2/„J") U Syzt„v), 

l<i<card(Gx) l<i<card(G„) 

where t £ Im-i is a bihomogeneous polynomial and the qi and q'^ are also bihomogeneous. Since 
^2 7^ and g^^^ £ k[xo, . . . ,XnJ, qi must be in {yo, - ■ ■ ,yny)- According to Proposition 14.31 

'^i,qig['^^ e Im-i- By a similar argument, Vi,?,-^!''^ £ Im-i- Finally, / £ Im-i- 

We just proved that Im-i ■ fm — Im-i U {Gx) U {Gy). Thus, G^; U U G,n-i is a Grobner basis of 

^m — 1 ■ fm • n 

CoroUarv 14. II shows that, when a bilinear system is bi-regular, it is possible to find a Grobner basis of 
Im-i ■ fm (which yields the monomials t such that the row {t, /,„) reduces to zero) as soon as we know 
the three Grobner bases Gx, Gy, and Gm-i- In fact, we only need Gx and Gy since the reductions to zero 
corresponding to Gm~i are eliminated by the usual F5 criterion. Fortunately, we can obtain Gx and Gy 
just by performing linear algebra over the maximal minors of a matrix (Theorem l3.2l l. 

We now present the main result of this section. If we suppose that Conjecture 14. II is true, then the 
following Theorem shows that generic bilinear systems are bi-regular 
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Theorem 4.1. Let m, n^;, rij, € N such that m < + riy. The set of bi-regular sequences (/i, . . . , /„) 
contains a nonempty Zariski-open set. Moreover, if {fi, ■ ■ ■ , fm) is a bi-regular sequence, then there are 
no reductions to zero with the extended criterion. 

Proof. Let G„i be a minimal Grobner basis of Im-i : fm- The reductions to zero which are 

not detected by the usual F5 criterion are exactly those such that t E LM(G„i) and t ^ LM(/„i_i). 
We showed that there exists a nonempty Zariski-open subset Oi of B£{nx, Uy) such that if /„, e Oi, 
then t e LM (/,„_! : /,„ n k[ ])ort e LM(/„_i : /„ n k[yo, . . . ,2/„J) (Proposition |42J. 

If we suppose that the conjecture 14.11 is true, then there exists a nonempty Zariski-open subset O2 of 
BC{nx, Uy) such that if /,„ e O2, I,n-i ■ fm n k[xo, • . • , XnJ (resp. : n k[yo, . . . , 2/„J) is 

spanned by the maximal minors of jaCx(i^m-i) (resp. jaCy(Fm_i)). Thus, by Theorem 13. 21 there exists a 
nonempty Zariski-open subset O3 of rij,) such that if /m e O3, LM(/„i_i : fmr\k[xQ, . . . , Xn^]) = 

MonomialsJ^j_„^_2('^a + 1)) (I'^sp- LM(/,„_i : /,„ n fc[?/o, • • • ,ynj) = Monomials^_„^_2(rix + I)))- 
Suppose that /,„ e Oi n O2 n O3 (which is a nonempty Zariski-open subset) and that {t, fm) is a reduction 
to zero such that t ^ LM(/,„_i). Then 

t e (Monomials^_„^_2(ny + 1)) 
or 

t G (Monomialsf„_„^_2(nx + 1)). 

By Lemma IX2I t is a leading monomial of a linear combination of the maximal minors of jac^{F„i-i) (or 
jaCy (F,„_i)). Consequently, the reduction to zero {t, fm) is detected by the extended F5 criterion. □ 

Remark 4.3. Thanks to the analysis of Algorithm 13.21 we know exactly which reductions to zero can be 
avoided during the computation of a Grobner basis of a bilinear system. If a bilinear system is bi-regular, 
then the Algorithm \3.2\f inds all reductions to zero. Indeed, this algorithm detects reductions to zero coming 
from linear combinations of maximal minors of the matrices ]ac^{Fi) and jaCy {Fi). According to Theorem 
\4.1\ there are no other reductions to zero for bi-regular systems. 

5 Hilbert bi-series of bilinear systems 

An important tool to describe ideals spanned by bilinear equations is the so-called Hilbert series. In the 
homogeneous case, complexity results for F^ were obtained with this tool (see e.g. |5 1). In this section, we 
provide an explicit form of the Hilbert bi-series - a bihomogeneous analog of the Hilbert series - for ideals 
spanned by generic bilinear systems. To find this bi-series, we use the combinatorics of the syzygy module 
of bi-regular systems. With this tool, we will be able to do a complexity analysis of a special version of the 
^5 which will be presented in the next section. 

We say that an ideal is bihomogeneous if there exists a set of bihomogeneous generators. The following 
notation will be used throughout this paper: the vector space of bihomogeneous polynomials of bi-degree 
{a, P) will be denoted by Ra,i3. If / is a bihomogeneous ideal, then la.p will denote the vector space 

/ni?a,/3. 

Definition 5.1 ( 04OI 1361 ). Let I be a bihomogeneous ideal of R. The Hilbert bi-series is defined by 

HS/(ti,t2)= J2 <^MRa,p/Ic.,l3)t'^4- 

Remark 5.1. The usual univariate Hilbert series for homogeneous ideals can easily be deduced from the 
Hilbert bi-series by putting ti = ^2 (see Ii36]l j. 

We can now present the main result of this section: an explicit form of the bi-series for bi-regular 
bilinear systems. 
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Theorem 5.1. Let /i ,...,/„ e i? fee a bi-regular bilinear sequence, with m < Ux + riy. Then 
where 

^^-in.+l)^^ _ + _ _ _ ^^-j, ^„^+i_fe^/+„,^_fe^j _ 

We decompose the proof of this theorem into a sequence of lemmas. 

If / is an ideal of R and / is a polynomial, we denote by / the equivalence class of f inR/I and 

annn/j{f) = {veR/I ■.vf = 0}, 

annii/i{f )a,/3 = {v e i?// of bi-degree {a,P) : vf = 0}. 
If 7 is a bihomogeneous ideal and / is a bihomogeneous polynomial, we use the following notation: 

Gij{ti,t2)= ^ dim(annfl/j(/)„,/3)fft^. 



Lemma 5.1. Let /i ,...,/„ e i? foe bihomogeneous polynomials, with 1 < m < + riy. Let {di ,d2) be 
the bi-degree of fm- Then 

Proof. We have the following exact sequence: 

^ ann^/,^_,(/) ^ R/Im-i ^ R/Im-i ^ R/Im ^ 0. 

where Lpi and (^3 are the canonical inclusions, and is the multiplication by /„. 

From this exact sequence of ideals, we can deduce an exact sequence of vector spaces: 



^ (ann«/,_^ ^ ( ) ^ ( ) ^ (f) ^ 0. 

Thus the alternate sum of the dimensions of vector spaces of an exact sequence is 0: 

dim((ann^/j^^_^(/))„,^) - dim (^{j^^ ^ ^ + 
dimff^^) Vdimfff) \=Q. 

\ \ -''"-1 / a+di,/3+d2 / \\'^J a+<ii,/3+d2 / 

By multiplying this relation by t^tj ^^id by summing over (a, /3), we obtain the claimed recurrence: 
HS/„(ti,t2) = {l-t'lHi^)BS,i^_,+t'iHi^Gi^_,j{hM)- 

□ 

Lemma 5.2. Let /i, . . . , S R be a bi-regular bilinear sequence, with m < Ux -\- riy. Then, for all 
2<i<m, 

where 

{0 ifi < Uy 
(l-t)"x + i 2^1<j<ny + l (i_t)"x+"a-i+2 
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Proof. Saying that v € ann^//.^ {fi) is equivalent to saying that the row with signature (LM(w), fi) is not 
detected by the classical criterion. According to Theorem |4.1| if the system is bi-regular, the reductions 
to zero corresponding to non-trivial syzygies are exactly: 



U {{t, fi) : t e Monomialsf^„^_2K + 1)} IJ f^)'-*^ Monomialsf_„^_2K + !)}■ 

By Proposition|43] we know that if P e k[xo, . . . , a;„^] n {li-i : fi) (resp. k[yo, . . . , ?/„J n {li-i : fi)), 
then UjP G (resp. xjP G li-i). Thus Gi-^-^j^ (^1,^2) is the generating bi-series of the monomials 
of k[xQ, . . . , Xn^] which are a multiple of a monomial of degree Uy + 1 in xq, . . . , Xi_„ _2 and of the 
monomials of k[yo, . . . , which are a multiple of a monomial of degree + I in yo, . . . , yi-n^-2- 
Denote by gi^~^\t) (resp. gy ^\t)) the generating series of the monomials of k[xQ, . . . ,Xn^] (resp. 
k[yo, . . . , yuy]) which are a multiple of a monomial of degree Uy + 1 (resp. + 1) in a;o, . . . , Xi-ny-2 
(resp. 2/0, ... , yi-n^-2). Then we have 

Gu_,j,{h,h)^gt'Hh)+gl;-^\t2). 

Next we use combinatorial techniques to give an explicit form of gx (t) and gy^'^'^ (t). Let c{t) denote 
the generating series of the monomials of k[xi^ny-i, ■ ■ ■ , Xn^]'- 

rix + Uy — i + j + 1 



j j~ (\ — fyix+''iy—i + 2' 

Let Bj denote the number of monomials of k[xo, . . . , Xi-ny-2] of degree j. Then 

(l„^)i+n,+2 = ^(t) + Bic{t) + ■■■+ B,,^c{t) + gt'Kt). 
Since Bj — (* ""j^^^"*), we can conclude: 

if i < n 



(l_f)nx + l Z^l<j<n„ + 1 (l_t)"a:+na-i + 2 



Proof of Theorem [5. 1 1 Since the polynomials are bilinear, by Lemma lSTl we have 

HS/.(ti,t2) = {l-tlt2)mu_,+ht2Gu_,J^{h,t2)- 

Lemma ISTSl gives the value of Gi^_^,f^ (^1,^2)- To initiate the recurrence, we need 
Then we can obtain the claimed form of the bi-series by solving the recurrence: 



(1 + -t2)"« + l 



N.,{tiM) = (1 - ht2y + tMi - ht2yGi^j^^,{h,t2) 

3=0 



□ 



□ 
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6 Towards complexity results 



6.1 A multihomogeneous F5 Algorithm 

We now describe how it is possible to use the multihomogeneous structure of the matrices arising in the 
Matrix F5 Algorithm to speed-up the computation of a Grobner basis. In order to have simple notations, 
the description is made in the context of bihomogeneous systems, but it can be easily transposed in the 
context of multihomogeneous systems. 

Let /i ,...,/,„ be a sequence of bihomogeneous polynomials. Then consider the matrices Aid in degree 
d appearing during the Matrix F5 Algorithm. One can remark that each row represents a bihomogeneous 
polynomial. Let (^1,^2) be the bi-degree of one row of this matrix. Then the only non-zero coefficients on 
this row are in columns which represent a monomial of bi-degree {di, ^2). Then a possible strategy to use 
the bihomogeneous structure is the following: 

• For each couple {di , ^2) such that di+d2 — d, construct the matrix Md^,d2 ■ The rows of this matrix 
represent the polynomials of Md of bi-degree {di, ^2) and the columns represent the monomials of 

• Compute the row echelon forms of the matrices Md^ .d2 ■ This gives bases of Id^ ,^2 ■ 

• The union of the bases gives a basis of /c; since ^ ®dl+d2=d^dl.d2■ 
T\\\s, way, instead of computing the row echelon form of a big matrix, we can decompose the problem 

and compute independently the row echelon forms of smaller matrices. This strategy can be extended to 
multihomogeneous systems. 

In Table [T] the execution time and the memory usage of this multihomogeneous variant of i^s are 
compared to the classical homogeneous Matrix Fc, Algorithm for computing a Z3-Gr6bner basis for random 
bihomogeneous systems (for the grevlex ordering). Both implementations are made in Magma2 . 15-7. 
The experimental results have been obtained with a Xeon processor 2.50GHz cores and 20 GB of RAM. We 
are aware that we should compare efficient implementations of these two algorithms to have a more precise 
evaluation of the speed-up we can expect for practical applications. However, these experiments give a 
first estimation of that speed-up. Furthermore, we can also expect to save a lot of memory by decomposing 
the Macaulay matrix into smaller matrices. This is crucial for practical applications, since untractability is 
often due to the lack of memory. 



6.2 A theoretical complexity analysis in the bilinear case 



In this section, we provide a theoretical explanation of the speed-up observed when using the bihomoge- 
neous structure of bilinear systems. To estimate the complexity of the Matrix F5 Algorithm, we consider 
that the cost is dominated by the cost of the reductions of the matrices with the highest degree. By using 
the new criterion described in Section 13.41 all the matrices appearing during the computations have full 
rank for generic inputs (these ranks are the dimensions of the fc-vector spaces /di.dg)- We consider that the 
complexity of reducing a ?- x c matrix with Gauss elimination is 0{r'^c). Thus the complexity of computing 
a Z?-Gr6bner basis with the usual Matrix F5 Algorithm and the extended criterion for a bilinear system of 
m equations over k[xo, . . . , a;„^ , yo, • ■ • , 2/« J is 



'- horn 





When using the multihomogeneous structure, the complexity becomes: 
/ 



multihorn 



= C2 



(dim(i?d,,dj- [tft^^]HS(ti,i2)) dim(i?d,, 

I di+d2 = D 

\ 1 < di, d2 < D - 1 
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Multihomogeneous 


Homogeneous 








m 


bidegree 


U 


time 


memory 


time 


memory 




3 


4 


7 


(1,1) 


6 


16.9s 


30MB 


265.7s 


280MB 


16 


3 


4 


7 


(1,1) 


7 


105s 


92MB 


2018s 


1317MB 


19 


4 


4 


8 


(1, 1) 


7 


582s 


275MB 


13670s 


4210MB 


23 


5 


4 


9 


(1, 1) 


7 


3343s 


957MB 


66371s 


12008MB 


20 


5 


5 


10 


(1, 1) 


6 


645s 


435MB 


10735s 


4330MB 


17 


2 


2 


4 


(1,2) 


10 


11.4s 


19MB 


397s 


299MB 


35 


2 


2 


4 


(1,2) 


8 


1.7s 


10MB 


16s 


52MB 


9 


3 


3 


6 


(1,2) 


8 


67s 


80MB 


1146s 


983MB 


17 


4 


4 


8 


(1, 2) 


8 


2222s 


1031MB 


40830s 


12319MB 


63 


2 


2 


4 


(2,2) 


11 


29s 


27MB 


899s 


553MB 


31 


3 


3 


6 


(2,2) 


8 


27s 


47MB 


277s 


452MB 


10 


3 


3 


6 


(2,2) 


9 


152s 


154MB 


2380s 


1939MB 


16 


3 


4 


7 


(2,2) 


9 


1034s 


505MB 


18540s 


7658MB 


18 


4 


4 


8 


(2,2) 


8 


690s 


385MB 


7260s 


4811MB 


11 


4 


4 


8 


(2,2) 


9 


6355s 


2216MB 




>20000MB 





Table 1 : Execution time and memory usage of the multihomogeneous variant of F5 



where dim(i?rfj ,i2) — ( ^^"^) ( ^j^"")- Thus the theoretical speed-up that we expect is: 



speedupth = CzF{nx, ny, m, D) 



where C3 = ^ is a constant and 



C2 



( 



F(nx,ny,m,D) 



(dim{Rd,A,) - [i?H^^]HS(ti,t2))'dim(i?d,,rf, 

1 di + d2 = D 

\ I<di,d2<-D-1 



Now let us compare this theoretical speed-up with the one observed in practice. 





Uy 


m 


D 


experimental 
speed — up 


F{nx,ny,m,D) 


3 


4 


7 


6 


16 


29 


3 


4 


7 


7 


19 


34 


4 


4 


8 


7 


23 


34 


5 


4 


9 


7 


20 


32 


5 


5 


10 


6 


17 


27 



We can see in this table that, in practice, experimental results match the theoretical complexity: 

speedup w 0.6F(nx, Uy, m, D). 

6.3 Structure of generic affine bilinear systems 

In this section, we show that generic affine bilinear systems have a particular structure: they are regular 
(Definition l2.6l l. Consequently, the usual F5 criterion removes all reductions to zero. 
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Proposition 6.1. Let S be the set of affine bilinear systems over k[xi, . . . , Xn^ , J/i, • ■ • , Vny] with m < 
nx + ny equations. Then the subset 

{(/i, • • ■ ,/m) e : {fi,...,fm) is a regular sequence} 

contains a Zariski nonempty open subset of S. 

Proof. Let (/i, . . . , /„) be a generic affine bilinear system. Assume that it is not regular Then for some 
i, there exists g £ R such that g ^ li^i and gfi S li-i- Denote by g'' the bi-homogenization of g. 
Then g^ G {f^, . . . , : //*. (/{', . • . , /^) is a generic bilinear system, hence it is bi-regular (Theorem 
14.1b . Thus g'^ G k[xo, . . . , XnJ\ or g^ S fc[yo: • • ■ : Vny]- Let us suppose that g^ G k[xQ, . . . , (the 
proof is similar if g'^ G k[yo, . . . , j/„j,]). Therefore Unyg^ G . . . , //Li) when the system is bi-regular 
(Proposition 14. 3l l. By putting a;„^ = 1 and ?/„ = 1, we see that in this case, g G li-i, which yields a 
contradiction. This shows that generic affine bilinear systems are regular. 

□ 



6.4 Degree of regularity of affine bilinear systems 

In this part, m, n^. and Uy are three integers such that m — rix + riy. We consider a system of bilin- 
ear polynomials F — {fi, . . . , fm) G k[xQ, . . . , Xn^ , yo, • ■ • , 2/™ ]™- denotes the deshomogenization 
morphism: 

fc[xo,...,x„^,yo,---,2/nJ — > k[xo,...,Xn^-i,yo,...,Vny-l\ 
f{xo,..-,Xn^,yo,.-.,y„y) ' — > /(a;o, . . . , x„^_i, 1, yo, • • ■ , yn^-i, 1) 

Also, / stands for the ideal (/i, . . . , /,„) and denotes the ideal (t?(/i), . . . , i9(/m)). In the follow- 
ing, we suppose without loss of generality that < riy. We also assume in this part of the paper that the 
characteristic of fc is (although the results remain true when the characteristic is large enough). 

The goal of this section is to give an upper bound on the so-called degree of regularity of an ideal / 
generated by a generic affine bilinear system with m equations and m variables. The degree of regularity is 
a crucial indicator of the complexity of Grobner bases algorithms: for 0-dimensional ideals, it is the lowest 
integer dreg such that all monomials of degree dreg are in LM(/) (see jSj). As a consequence, the degrees 
of all polynomials occurring in the Fc, algorithm are lower than dreg + 1. In the following, -< still denotes 
the grevlex ordering. 

Lemma 6.1. If the system F is generic, then there exists polynomials go, . . . ,gn^-i G k[yo, . . . , y„^_i] 
such that 

Vj G {0,. . - l},Xj - gj{yo, . . . ,y„^-i) G ■&{!). 
Proof. We consider the m x rix matrix A = iac^{'d{F)) and the vector 

B - (z9(/i)(0, . . . , 0, yo, . ■ . , Vny-i) ■ • ■ ^?(/™)(0, . . . , 0, yo, . . . , y„,-i)) . 





/ Xo \ 






Thus A ■ 




+ B = 





We denote by {A^*)} all the rix X Tlx sub-matrices of A. 
Let (ao, . . . , a„ _i) G V ar{{MaxM\nors{ii}Qac^{F))))) be an element of the variety. Let yl^ (resp. i?^) 
denote the matrix A (resp. B) where y^ has been substituted by ai for all i. Since is 0-dimensional, 
the affine linear system 

A„ • [ ... j + B„ ^ 

has a unique solution. Therefore, the matrix Aa is of full rank. Consequently, there exists an invertible 
n,, X n„ sub-matrix of A„. 
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Since k is infinite, we can suppose without loss of generality that, if the system is generic, then for all 
a, the matrix Aa'' obtained by considering the Uy first columns of Aa is invertible (if Aa^ is not invertible, 
just replace the original bilinear system by an equivalent system where each new equation is a generic 
linear combination of the original equations). Thus det{Aa^) ^ 0. 

According to Lemma[33]and|B31 (MaxMinors(i?(jac^(i^)))) = . . . , i?(/™))nfc[2/o, • • • ,?/n„-i]. 

Thus det(A'^^') (i.e. the matrix of the riy first columns of A) does not vanish on any elements of the va- 
riety of Therefore, the Nullstellensatz says that det(A*-^^) is invertible in fc[yo, ■ ■ ■ ^Vn H 
fc[yo, • ■ • 1 Vny-i])- Let h denotes its inverse. We know from Cramer's rule that there exists polynomials 
gj ek[yo,..., yny-i] such that 

Xj det(A(i)) - .g,(yo, • ■ • , Vny-i) e 
Multiplying this relation by h, we obtain: 

Theorem 6.1. If the system F is generic, then the degree of regularity of'd{I) is upper bounded by 

dreg < min(nj, + l,ny + 1). 

Proof. We supposed that < Uy, so we want to prove that dreg = rix + l- Let t = YYj=o^ ^'j^ Ofc^o ^ Vk^ 
be a monomial of degree rix + I- According to Lemma l6n 

n^-l riy-1 

t- n 9jiyo,-.-,yny-ir n yt^^i^)- 

j=0 k=0 

Now consider the normal form with respect to the ideal (MaxMinors(i9(jaCx(-F)))). Then 

j=0 k=0 

Since all monomials of degree Ux + 1 are in LM((MaxMinors(??(jaCx(F))))) (Lemma [3.2l) . 

deg(NF( n 9j{yo, . . . , Vny-ir n 2^^)) < + 1- 

j=0 k=0 

This implies that 

LM(i-NF(n 5,(2/0,..., J/„„-ir n yk'))^t- 

j=0 k=0 

Therefore, for each monomial t of degree rix + l,t E LM(i?(/)). This means that dreg < + 1- □ 

Remark 6.1. This bound on the degree of regularity should be compared with the degree of regularity of 
a generic quadratic system with m equations and m variables. The Macaulay bound ( see Ii26\l } says that 
the degree of regularity of such systems is m + 1. Since Grobner bases algorithms are exponential in the 
value, it means that ajfine bilinear systems are a lot easier to solve than generic ajfine quadratic systems. 
Moreover, the inequality dreg < min(nj; + 1, rij, + 1) is sharp: experimentally, it is an equality for random 
bilinear systems. 

The following Corollary is a consequence of Theorem 16. II 

Corollary 6.1. The arithmetic complexity of computing a Grobner basis of a generic bilinear system 
/i , . . . , fn^+uy e fc [xo , . . . , _ 1 , 2/0 , . . . , yn^ - 1 ] With the F5 Algorithm is upper bounded by 

fix + riy + min(n2^ + l,ny + 1)\ "\ 
min(nj; + l,ny + 1) J J ' 

where 2 < uj < 3 is the linear algebra constant. 
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Proof. According to fS\, the complexity of the computation of the Grobner basis of a 0-dimensional ideal 
is upper bounded by 




where n is the number of variables and dreg denotes the degree of regularity. In the case of a generic affine 
bilinear system in k[xQ, . . . , Xn^-i,yo, ■ ■ ■ ,yn -i], n = Ux + riy and dreg < min(nj; + l, n^ + l) (Theorem 

lO) . ' " □ 



7 Perspectives and conclusion 

In this paper, we analyzed the structure of ideals generated by generic bilinear equations. We proposed an 
explicit description of their syzygy module. With this analysis, we were able to propose an extension of 
the F5 criterion dedicated to bilinear systems. Furthermore, an explicit formula for the Hilbert bi-series is 
deduced from the combinatorics of the syzygy module. With this tool, we made a complexity analysis of a 
multihomogeneous variant of the Algorithm. 

We also analyzed the complexity of computing Grobner bases of affine bilinear systems. We showed 
that generic affine bilinear systems are regular, and we proposed an upper bound for the degree of regularity 
of those systems. 

Interestingly, properties of the ideals generated by the maximal minors of the jacobian matrices are 
especially important. In particular, a Grobner basis (for the grevlex ordering) of such an ideal is a linear 
combination of the generators. In the affine case, this ideal permits to eliminate variables. 

The next step of this work would be to generalize the results to more general multihomogeneous sys- 
tems. For the time being, it is not clear how the results can be extended. In particular, it would be interesting 
to understand the structure of the syzygy module of general multihomogeneous systems, and to have an 
explicit formula of their Hilbert series. Also, having sharp upper bounds on the degree of regularity of 
multihomogeneous systems would be important for practical applications. 
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A Bihomogeneous ideals 

In this part, we use notations similar to those used in Section]?) 

• BT-l{nx, riy) the fc- vector space of biUnear polynomials mk[xQ, . . . , ,yo, . . . , yuy]', 

• X C k[xQ, . . . , , yo, ■ • ■ , J/n„] (resp. Y) is the ideal ( ) (resp. (yo, ■ ■ • ,2/«J); 

• An ideal is called bihomogeneous if there exists a set of bihomogeneous generators. In particular, 
ideals spanned by bilinear polynomials are bihomogeneous. 

• Ji denotes the saturated ideal U : (X n 

• Given a polynomial sequence (/i , . . . , /„), we denote by Syztriv the module of trivial syzygies, i.e. 
the set of all syzygies (si, . . . , Sm) such that VI < i < m, s, e {fi, . . . , fi-i,fi+i, /„); 

• A primary ideal P C i? is called admissible if {xq, ■ ■ . , Xn^) (t \/P and (yo, ■ • ■ , J/n„) VP', 

• Let £' be a fc-vector space such that dim(i?) < oo. We say that a property V is generic if it is satisfied 
on a nonempty open subset of E (for the Zariski topology), i.e. 3h E k[ai, . . . , adim(£;)], h 0, 
such that 

P does not hold on (ai, . . . , adim(_E)) Hai, . . . ,adim(£;)) = 0. 
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Proposition A.l ( 11361 ). Let I be an ideal of R. The two following assertions are equivalent: 

• / is bihomogeneous. 

• For all h Cz I, every bihomogeneous component ofh is in I. 

Lemma A.l ( II36I ). Let /i, . . . , /,„ £ Rbe polynomials, and /,„ = C\Pi be a minimal primary decompo- 
sition of Im ond let Adm be the set of the admissible ideals of the decomposition. Then Jm = (^PeAdmP- 

Proposition A.2. let /i, . . . , ^ R be polynomials with m < n^. + Uy, and Ass{Ii_i) be the set of 
prime ideals associated to li-i. The following assertions are equivalent: 

L V2 < i < m, fi is not a divisor of Q in Rj Ji^\. 

2. \/2 <i < m, [fi E P,P E Ass{Ii-i)) P is non-admissible. 

Proof. It is a straightforward consequence of Lemma lATTI □ 

Remark A.l. All results in this section can be generalized to multihomogeneous systems. Since we focus 
on bilinear systems in this paper, we describe them in this more restrictive context. 

Lemma A.2. Let P be an admissible prime ideal of R. The set of bilinear polynomials f (z R such that 
f ^ P contains a Zariski nonempty open set. 

Proof. Let / be the generic bihnear polynomial 

in k{{aj^k}vi<j<n^fi<k<ny)[xo, . . . ,Xn^,yo, . . . , yn^]- Since P is admissible, there exists xj^yko such that 
^joVko ^ P (this shows the non-emptiness). Let -< be an admissible order. Then consider the normal form 
for this order 

NFp(/) = ^t(oo,o---,a„,,„Jt. 

t monomial 

where the /if 's are polynomials. Thus, if a polynomial / e i? is in P, then its coefficients are in the variety 
of the polynomial system Vi, /it(ao,o, ■ ■ ■ , ^n^cny ) = 0. □ 

Theorem A.l. Let m,nx,ny € N such that m < rix + Uy. Then the set of bilinear systems /i, . . . , fm 
such that V2 < I < m, fi is not a divisor o/O in R/ Ji^i contains a Zariski nonempty open subset. 

Proof. We prove the Theorem by recurrence on m. Suppose V2 < i < m — 1, is not a divisor of 
in Rj Ji-\. We prove that the set of bilinear polynomials / such that / is not a divisor of in R/ Jm-i 
contains a nonempty Zariski open subset. According to Lemma lA.21 for each admissible prime ideal 
P e Ass{Im-i), the set Op — {f ^ P} contains a nonempty Zariski open subset. Thus f]p Op contains 
a nonempty Zariski subset (since the intersection of a finite number of nonempty Zariski open subsets is a 
nonempty Zariski open subset). Therefore, the set of bilinear polynomials / which are not divisor of in 
R/ Jm-i (this set is exactly f]p Op) contains a Zariski nonempty open subset. □ 

Proposition A.3. Let m < rix + riy and /i , . . . , be bilinear polynomials such that V2 < z < m, /i is 
not a divisor ofO in R/ Ji-i. Then VI < z < to, the ideal Ji is equidimensional and its co-dimension is i. 

Proof. We prove the Proposition by recurrence on to. 

• Ji = /i is equidimensional and codim(/i) = 1; 

• Suppose that Ji_i is equidimensional of co-dimension i — 1. Then Ji — ( Ji_i + fi) : {X D Y)°°. 
fi is not divisor of in Ji_i (Theorem I A. lb . thus Ji_i + fi is equidimensional of co-dimension 
i. Next, the saturation does not change the dimension of any primary component of a minimal 
primary decomposition of Ji_i + fi (the saturation only removes some components). Therefore, Jj 
is equidimensional and its co-dimension is i. 

□ 
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B Ideals generated by generic affine bilinear systems 

Let fc be a field of characteristic 0, m = + Uy, and a be the set 

o — \ : 1 < « < "^, < j < na;, < k < Uy}. 
We consider generic polynomials /i, . . . , /„ in /c(a)[a;o, . . . , , yo, • ■ • , ynj: 

and we denote by / C A:(a)[a:;o; . . . , ,yo, . . . ,yn ] the ideal they generate. In the sequel, denotes the 
deshomogeneization morphism: 

k[xo,...,Xn,,yo,.-.,yny] — > k[xa,...,Xn^-i,yo,...,yny-i] 
f{xo,.--,Xn^,yo,---,y„y) ' — > /(a;o, . . . , x„^_i, 1, j/o, ■ • ■ , yn„-i, 1) 

For a £ fc™("x+n»+2)^ stands for the specialization: 

: k{a)[xo,...,Xn^,yo,...,ynJ A;[a;o, . . . , , yo, ■ • • , 2/nJ 

f{a){xo,...,Xn^,yo,...,yny) f{a){xo,...,Xn^,yo,...,yny) 

Also Var{ipg,{I)) C P""= x P"* (resp. l/ar(?? o ipsi{I)) C denotes the variety of (pa{I) (resp. 

Lemma B.l. There exists a nonempty Zariski open set Oi such that if a. G Oi, then for all [ao, . . . , , 
/3o, • • ■ , Priy) G Var{(pii{I)), a.n^ =/= one/ =/= 0. TTz/i implies that the application 

Var{i9 o ipail)) — ^ T/ar(^a(/)) 

(ao, ■ • ■ ,an^-i,^o, ■ • ■ ,/3n„-i) I — > (ao, • ■ • , a„^-i, 1, /3o, • ■ • 1) 

/i fl bijection. 

Proof. See EOl page 75 1 ] . □ 

Lemma B.l. There exists a nonempty Zariski open set O2, such that if a. £ O2, then the ideal 1} o ipa{I) is 
radical. 

Proof. Denote by F the polynomial family (/i, . . . , /„). Let J C k[a] be the ideal (/ + (det(jac(i^)))) n 
k[a] and ^ be its associated algebraic variety. By the Jacobian Criterion (see e.g. |12, Theorem 16.19]), 
if a does not belong to then o ipa{I) is radical. Thus, it is sufficient to prove that fc™("^+"B+2) ^ ^ 
is non-empty. 

To do that, we prove that for all a e i^m{n^+ny+2) ^ there exists (ei, . . . , e„i) such that the ideal (t? o 
fa.{fi) + El, . • . , i9 o ifaifm) + £m) IS radical. Denote by gi ^ -d o (paifi) for 1 < i < m and consider the 
mapping 5' 

xefc™^(gi(x),...,g™(a;))efc™. 

Suppose first that \E'(fe'") is not dense in fc™. Since 'ii{k™ ) is a constructible set, it is contained in a Zariski- 
closed subset of fc™ and there exists (£!,...,£„) such that the algebraic variety defined by gi — ei = 
• ■ ■ = 9m — £m = is empty. Since there exists a' such that gi — £i = ■& o (fa'ifi), we conclude that 
-do ipa'{I) = (1). This implies that a! ^ J . 

Suppose now that "^(k™^ is dense in fc™. By Sard's theorem fEV, Chap. 6, Theorem 2], there exists 
(£!,...,£,„) e fc™ which does not lie in the set of critical values of ^. This implies that at any point of the 
algebraic variety defined by gi — £1 = • • • = g,„ — £,„ = 0, i9 o (^a(det(jac(F))) does not vanish. Remark 
now that there exists a' such that gi — £,; = o Lp^' (fi). We conclude that a' e i^rn{n^+ny+2) ^ which 
ends the proof. □ 

Lemma B.3. There exists a nonempty Zariski open set O3, such that //a G O3, 



^(MaxMinors(i? o ^a(jaCy (F)))) = (i? o ^^(/i), • • • , o (^a(/m)) n fc[ 
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Proof. Let a be an element in O2 (as defined in Lemma IB. 21 ). Thus d o ipa_{I) is radical. Now let 

{vo, ■ ■ ■ , -i,wo, ■ ■ ■ , - 1 ) G Var{'d o If he an element of the variety. Then 



(1? o (^a(jaCy(F))^,=t,^ 



f ^0 \ 



Wny-l 

V 1 / 




This imphes that rank(i9 o i^a(jaCy(F))j;.=„. ) < riy + 1, and therefore 

(wo, • ■ • , Wn^-i) e l/ar((MaxMinors(z9 o (pa(jaCy(i^))))). 

Conversely, let (wo, . . . , Wn^-i) G yar((MaxMinors('!9o(pa(jaCy (F))))). Thus there exists a non trivial 
vector (wo, • • ■ , u^rij,) in the right kernel Ker(i? o (pa(jaCy (F))^;.^^.). This means that (wq, ■ • ■ , Vn^-i, 1, 
Wo, ... , w„j^) is in the variety of 



{vq, . . .,Vn^-i, l,wo, . . .,Wny) € Var{ip^ (jaCy(F)) 



From Lemma lBTl 7^ if the system is generic. Hence 

{vo,.-.,Vn,-i, ) e Var(§o ip^{I)). 

Finally, we have 

Var{{MaxMmors{d o (^a(jaCy (F))))) = l^ar((i? o (^a(/i), • • • , o y^al/m)) n fc[xo, 
and o ifa{I) is radical (Lemma [B.2l l. The Nullstellensatz concludes the proof. 
Corollary B.l. There exists a nonempty Zariski open set O4, such that if a. e O4, 



' "^^j^ — 1]) 



□ 



Proof. According to Lemma IBT2] and Lemma lBTl if a e Oi n O2, then deg('!? o (pa.{I)) = card{Var{'d o 
fsi{I)) = card(Far((/3a(-f)))- This value is the so-called multihomogeneous Bezout number of ipa.{I), i-e. 
the coefficient of z"^ i" i^i + 22)""=+"" (see e.g. ill), namely □ 

Remark B.l. Actually, by studying ideals spanned by maximal minors of matrices whose entries are linear 
form, it can be shown that, for a generic affine bilinear system, (MaxMinors(i? o (pa(jaCy (F)))) is radical 
(see Lemma 1331 ). Hence Lemma \B.3\ shows that, for generic affine bilinear systems, 

(MaxMinors(i9 o ^a(jaCy(F)))) = (1? o ^a(/i), o ^a(/m)) n k[ 

(MaxMinOrs(l? O ^a(jaCx(F)))) = {§ O (^a(/l), • • • , O Lp^{frn)) n fc[2/o, ■ • • , Vny-l]- 



32 



